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We study the mean field equation of motion for driven domain walls in random media. We discuss 
the two cases of an external constant as well as an oscillating driving force. Our main focus lies on 
the critical dynamics close to the depinning transition, which we study by analytical and numerical 
methods. We find power-law scaling for the velocity as well as the hysteresis loop area. 
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I. INTRODUCTION 

The interaction of elastic systems, like charge density 
waves [3, 0] , Wigner crystals jj] , domain walls 0-0] , dis- 
location lines or magnetic flux lines with disorder is a 
problem of great technological importance. 

More than two decades ago, in a seminal work, D.S. 
Fisher has studied the depinning of charge density waves 
(CDW) from randomly distributed pinning centres [H H| 
by an external constant (dc) field h. He showed that the 
depinning transition is a dynamical critical phenomenon 
where the velocity close to the depinning threshold h p 
plays the role of an order parameter exhibiting a power 
law behaviour v ~ (h — hn)^ . Within the mean- field ap- 
proach, developed in [l|, Q, the exponent 8 was found 
to be 8 — 3/2. In a subsequent work Q, the treat- 
ment was extended to incorporate fluctuation effects us- 
ing a functional renormalisation group (FRG) approach. 
The latter is necessary since the correlator of the random 
pinning forces develops a cusp singularity on sufficiently 
large length scales. This results in a modified mean field 
exponent 8 = 1, to which corrections in e = 4 — d were 
calculated in 

A closely related problem concerns the behaviour of 
driven interfaces in random environments, as for exam- 
ple domain walls in disordered ferroic systems [l(| [H| 
or interfaces between immiscible fluids that are pushed 
through porous media [TU, [l3| . Feigel'man has con- 
sidered the velocity corrections due to the disorder and 
estimated from that the depinning threshold. However, 
taken self-consistently, the theory does not yield a finite 
depinning threshold. Koplik and Levine [l4| used pertur- 
bation theory for the associated mean-field model, which 
they derived analogous to the procedure by Fisher for 
CDWs [l[ ■ They found that the interface either follows a 
solution which moves with constant velocity or it remains 
pinned. They were, however, not able to extend their 
findings to the problem of spatially extended interfaces, 
because their perturbativc approach lacks the necessary 
FRG analysis. In a subsequent work, Leschhorn studied 
the mean field theory for domain walls in a model which 
treats the disorder in a simplified manner (lo| . He consid- 
ered a discretised lattice system and allowed the random 
force field to take values out of three possibilities only: 
— 1, or +1. For his model, he also found pinned and 



sliding solutions and determined the velocity exponent as 
8 = 1, which is the same as for CDWs when the disorder 
force has discontinuous jumps . Vannimenus and Der- 
rida (l6j simplified the Leschhorn model even further and 
were able to derive an exact solution. The basic simpli- 
fication of their model concerns the assumption of unit 
moves. This means, that per unit time step a segment of 
the interface either remains at rest, if the total force is 
smaller or equal to zero. Otherwise, it moves exactly one 
step forward, independent of the magnitude of the force. 
Though this assumption admits an exact solution, the 
restriction to unit moves entails a non-uniform periodic 
behaviour of the mean velocity close to the threshold. 
The time averaged velocity (over one period) has then a 
different exponent 8 = 1/2. 

In 1992, the depinning transition of interfaces in ran- 
dom media has been analysed within an FRG calcula- 
tion, starting from an expansion around a flat interface 
[171 fig . Similar to the charge density wave case, a cusp 
singularity of the force correlator develops giving rise to 
a finite depinning threshold as well as a non-trivial ex- 
ponent 8 = 1 — e/9, where d = 4 — e is the interface 
dimension. Later, this FRG calculation as been refined 
to include also two- loop terms [ill H^j ■ 

The FRG flow equations for dc driven interfaces have 
then been extended to the case of an oscillating (ac) driv- 
ing force [2l[. Combined with a scaling analysis it was 
possible to work out the main characteristics of the veloc- 
ity hysteresis loop v{h). In the limit of small frequencies, 
scaling behaviour has been found. The exponents of the 
remanent velocity at depinning as a function of frequency 
have been determined for all dimensions d < 4 and the 
results agree very well with the values obtained from a 
numerical study [2l| . The problem of ac-driven interfaces 
is also experimentally relevant. Of special interest is the 
ac susceptibility of ferroic systems [22], [23j] which gets a 
considerable contribution from the domain wall motion 
p4l [25| . A phenomenological understanding of different 
regimes has been reported in Ref. [26| . where the con- 
cept of waiting time distributions has been used. More- 
over, the perturbation theory for ac-driven interfaces in 
random environments has been examined [27l |. Further 
study of ac-driven elastic systems in disordered media 
has been devoted to vortex lattices [28| and structural 
defects in liquid crystals p9l |. 

The aim of this work is to study the properties of the 
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mean field theory for driven domain walls in random en- 
vironments close to the depinning threshold for dc- as 
well as for ac-driving forces. Especially for the regime 
of small driving frequencies the critical behaviour has 
not yet been investigated. We find scaling behaviour for 
small frequencies. We separately discuss the situation 
of smoothly correlated randomness, which is the usual 
starting point for any kind of FRG study of clastic sys- 
tems and a type of disorder the correlator of which has 
a cusp singularity at the o rig in, reflecting the fixed point 
solution of the FRG flow 0,^. 

The present article is organised in the following way. 
After a brief derivation of the mean field equation of 
motion in section [Til we focus on the special case of 
an adiabatic driving force in section IIIII We general- 
ize the arguments of Fisher Q to the problem of non- 
periodic systems, in order to find the scaling properties 
close to the depinning transition for adiabatic driving. 
Our analytical findings are supported by numerical anal- 
yses. In section [TV] we pay special attention to the small 
frequency scaling of several quantities characterising the 
velocity hyteresis loop. In the case of ac-drivings, ana- 
lytical treatments are difficult, so we extensively have to 
resort to numerical investigations. An outline about the 
numerical methods employed is provided in appendix iBl 



II. EQUATION OF MOTION 

The equation of motion of a driven interface in a d+ 1- 
dimcnsional inhomogeneous medium can be written as 



.10] 



1 dz 
7 dt 



T V 2 z + h(t) + u g(x,z), 



(1) 



where z(x,t) denotes the height of the interface profile. 
The three terms on the right hand side of equation ([T]) 
denote the curvature force, the external driving force 
h(t) = ho cos(tJoi) and the random force uog(x, z) aris- 
ing from the inhomogeneities of the medium. The coef- 
ficient 7 denotes the bare mobility. In the present paper 
we consider a mean field version of this equation. Mean 
field theory is suitable to describe systems above the up- 
per critical dimension. Below this critical dimension, the 
mean field analysis with a disorder correlator that has a 
cusp singularity at the origin (see below) yields the zcroth 
order for the critical exponents of the depinning transi- 
tion in an e-expansion (e = D c — D) [3j. In App. [A] 
we discuss the applicability of mean field theory in more 
detail. To obtain the mean field equation, we replace 
the local curvature term by some long range interaction 
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■£(z(t)-z(xj), (2) 



where y and x are nearest neighbour lattice sites, a is the 
lattice constant and 



z(t) 



d J2 z ( x ) ->■ L ~ d J d d x z ( x > *)■ ( 3 ) 
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in the continuum limit. This replacement results in an 
effective equation for z(x,t) 

1 dz(x,t) 2ci f—r-r , s\ , / x , , x x 
^ = ^r (z{t) - z{x,t))+h(t)+u g(z(x,t),x), 



(4) 

where x plays now the role of a label for z(t) and g(z). In 
the thermodynamic limit, L — > oo, z(t) does not fluctuate 
and hence (5]) can be solved for a given z(t) with the 
condition ([3]). For each x the solution depends on the 
random force field configuration {g(z)} for this particular 
value of x. The average over all positions, L~ d J d d x 
can finally be replaced by the average over the different 
random force configurations J Dg(z)P(g(z)), i.e. 



z(t) 



Dg(z)P(g(z))z(t) = (z(t)) . 



(5) 



For the random force we assume a Gaussian distribution 
with (g(z)} = 0, and 



(g{z)g{z')) = A e (z - z'), 



(0) 



where A^(0) = 1 and / A^(z)dz ~ £. Here, I denotes the 
correlation length of the force correlator, i.e. Ai(z) <C 1 
for \z\ ^> I. In this paper we are going to consider two 
different types of correlators. We distinguish between a 
correlator that is smooth and a correlator that shows a 
cusp singularity at the origin. 

In the following we measure z in units of I and t in 
units of £/"fUg, such that 

dz 

— =T((z) - z) + hcosojt + g(z), (7) 
at 

where T = 2dT $t / (uqcl 2 ) , uj = oJo£/(~fiio) and h = ho/uo. 
The function g(z) obeys the same relations as above with 
i = 1, and therefore we will skip the index I at A^(z) 
from now on. Thus, the theory depends on the three pa- 
rameters r, h and u>. The mean field equation of motion 
is similar to that for charge density waves considered by 
Fisher [3(| if g{z) is replaced by cos(z — j3) where /3 is a 
random phase. 

The physical picture of the mean-field equation of mo- 
tion is, that segments of the interface now behave like 
individual particles, each of which moves in a distinct 
configuration of the disorder. Every particle is coupled 
to the disorder averaged position, which in turn is deter- 
mined self-consistently. 



III. ZERO FREQUENCY LIMIT 

In this section, we consider a special case of the equa- 
tion of motion ((7J) , for which the driving force is constant 



3 




FIG. 1: Plot of the left and the right hand side of equation 
(|13|l for random forces with a smooth correlator. For small 
F there are several intersection points, for small values of V 
there is only one solution. 



in time 



dz 

at 



T((z) - z) + h + g{z). 



(8) 



At sufficiently large driving force h, the average particle 
position (z) will move with constant velocity v > 0. In 
this case, Eq. ((8]) can be written as 



d t z = Tvt + h-Tz + g(z) 
= g{z)-[T(z-vt)-h}. 



(9) 



In the following we will consider the case where the ve- 
locity is sufficiently small v -C h. The positions where 
d t z = follow from the intersection of g{z) with the 
straight line Y{z — vt) — h which moves to the right with 
velocity v (cf. Fig. [1]). For sufficiently small F and 
smooth g{z) there are in general In + 1 intersections 
which we denote by z\ < Z2 < z% < . . .. For z < Z2, z is 
driven towards z\ , for zi < z < z± it is driven towards zj, 
etc. If the particle starts with an arbitrary initial value, 
it will first develop towards the closest stable fixed point 
of ©, where the particle velocity is small. Let us as- 
sume this is Z\. The force free point z\ will then change 
according to z\(t) = vt + h/T + g(zi)/T. Eventually, 
the intersection point z\(t) merges with zi{£) and then 
disappears. In this case z(t) will grow sufficiently fast un- 
til it reaches z^{t) and the process repeats if we replace 
z n — > z n -2- Thus, the motion of the particle is jerky: 
periods of slow motion with velocity v are intermitted by 
fast periods where the particle is driven towards a new 
stable fixed point. Below, we will analyse this process in 
detail. 



A. The general picture 

A first overview results from considering some limiting 
cases. 

(i) For large but finite T we can apply perturbation 
theory. To lowest non-trivial order one obtains for the 
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h p (T ) (for cusped disorder, num.) 
Linear fit fcpfT- 1 ) = 0.96 ■ T" 1 for small T" 1 
/ip(r _1 ) (for smooth disorder, num.) 
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FIG. 2: Depinning threshold as a function of T~ in the case 
of a dc-drive, ui — 0. For the case of a cusp-correlator of the 
random forces (diamonds) the depinning threshold remains 
finite as long T _1 is finite. For a smooth correlator (crosses) 
the threshold vanishes for small T _1 as expected from pertur- 
bation theory. 



mean velocity (a short derivation of this expression can 
be found in appendix |C|) 



v = h + 



/ dter Tt ^{vt). 
Jo 



(10) 



The depinning threshold h p ^± for h ^ follows from 
taking the limit v — > ±0, h — > ±h p ± 



h p ± = - limr^A'rie). 

£->0 



(11) 



Thus, the force correlator has to have a cusp singularity 
to produce a finite threshold. If there is no cusp, pertur- 
bation theory in signals the absence of a depinning 
threshold. This argument applies however only to the 
region where perturbation theory is applicable, i.e. for 
r ^> 1. This pcrturbative result is in accordance with 
our numerical analysis, as is shown in Fig. [2j 

As has been mentioned in the introduction, a cusp sin- 
gularity in the correlator emerges as a fixed point solution 
of the functional renormalization group (FRG) flow in 
4 — e dimensions and describes the effective randomness 
on scales larger than the Larkin length. This leads to 
the existence of a depinning threshold in all dimensions 
d < 4. Of course, in the framework of the mean field 
approximation an FRG study is senseless and a correla- 
tor with a cusp singularity has to be included manually. 
Nevertheless, as we already see here, in many aspects 
the assumption of a correlator with a cusp gives different 
results compared to a smooth correlator. Incidentally, 
Narayan and Fisher [j| |3l| have used the mean field so- 
lution for cusped disorder to expand around in order to 
work out the critical behaviour of 4 — e-dimcnsional sys- 
tems. 

(ii) Finally, we consider the case f< 1. For T = the 
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FIG. 3: Plot of the left and the right hand side of equation 
(|13|) for random forces accociated to a scalloped potential, 
which shows a cusp singularity in the correlator. Random 
force realisations with a jump close to the origin yield more 
than one solution even for very small values of F. 



equation of motion ((SJ can be integrated 

f z dz' , 4 

To calculate the integral we assume that h > and h + 
g(0) > 0. Then for small z the lefthand side is positive 
and hence t as well, so the velocity is finite. However, 
since g{z) is unbounded there is a value z\ at which the 
denominator vanishes. The integral is then dominated by 
the integration in the vicinity of z\ Thus if z approaches 
z\ the time scale t diverges and the velocity vanishes, the 
particle is pinned at z = Z\ . The same argument works 
for h<0. 



depinning all particles are localised at the leftmost force 
free points. 

Let us now apply this argument to the case r > 1. 
For a smooth potential as depicted in Fig. [1] there is 
typically only one solution z[. For this single solution, 
g(z[) can be positive or negative with equal probability. 
Thus, a pinned solution obeying h = — (g(z[)) does not 
exist apart from the case h = 0. Hence, in this case the 
interface is never pinned, in agreement with our result 
from perturbation theory. The situation is different in 
the case when the random force exhibits infinite slopes 
as is shown in Fig. [3J Then, due to the discontinuities 
there are in general several solutions z.; for any value of 
r from which the leftmost ones dominate the behavior in 
the neighborhood of the depinning threshold. Of course, 
the larger the value of T the smaller is the fraction of 
disorder realisations which allow for more than one force 
free solution. Thus, for large T, the depinning threshold 
is diminished but finite. The presence of infinite slopes 
is a special feature of disorder forces, the correlator of 
which has a cusp singularity at the origin. A detailed 
analysis of the cusped disorder, as is sketched in Fig. [3]is 
presented in App. [D] where we discuss how such a class 
of disorder forces can be realised and derive the correlator 
explicitly. 

To test these predictions, we have solved equation ([8]) 
numerically. The depinning threshold h p is plotted in 
Fig. H as a function of T. It is clearly seen that the 
threshold increases with T -1 , it vanishes for r > T c for 
smooth random force correlations. For cusp correlations 
r c — > oo. These findings support the results from per- 
turbation theory. 



B. Static solution 

One special class of solutions to the equation of mo- 
tion ([SJ are the static solutions z s with d t z s = 0. Here, 
we are going to analyse under which circumstances such 
solutions can exist |32| . 

From the equation of motion Eq. (|SJ it is clear that 

T(z s -(z))-h = g(z s ) (13) 

must be obeyed, i.e. the system has to be located at 
force- free positions. Besides Eq. (jTB")) . one has to take 
into account that the self-consistency condition 

h=-(g(z s )), (14) 

which follows from averaging (|13l) holds. The maximal 
value on the righthand side of ([T4")) is realised, if z s = z\ . 
Thus, 

h p = - (g( Zl )) (15) 

is a critical field strength, above which no static solutions 
are possible. Conversely, we can conclude that close to 



C. Scaling behaviour above depinning 

Now, we consider the behaviour slightly above the de- 
pinning threshold h > h p , when (z) = vt but »<1, Our 
goal is to work out the scaling exponent /3 for the sliding 
velocity v, which we anticipate to vanish as a power law 

v~ (h-h p f. (16) 

To this aim, we solve the equation of motion in an ap- 
proximate manner. As the velocity of the interface is 
small, v <ti 1, we can also expect that dtz <^ 1 for most 
of the time. Thus, z(t) follows essentially from the van- 
ishing of the righthand side of ([5]), which means that z(t) 
stays close to the leftmost fixed point zi(t). Since the 
disorder averaged position (z) is in motion, we have to 
keep in mind, that the root of the straight line in Figs. 
[Hand [3] is now moving relative to g(z). The intersection 
point zi(t) satisfies the relation 

«i(t) = ^ + r- 1 (/ l + .g(z 1 (t))). (17) 

Without loss of generality, we restrict ourselves to v > 
0, so zi(£) moves now to the right. In this part of the 
motion, zi changes slowly (of order v). Eventually, z\ 
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merges with z<i. Let us assume that this happens at t = 0. 
For further reference we denote 



'-0 



L (0-)=z a (0-). 



(18) 



For t > 0, these two solutions disappear and the intersec- 
tion point 23(0—) becomes the new leftmost intersection 
point, i.e. 23(0—) — > zi(0+), so effectively z\ jumps in- 
stantaneously. Thus, at t > 0, the position z(t) is not 
any more close to a force free position and therefore it 
moves faster to approach the new intersection point z\ (t). 
The idea is now, that the mean velocity v is mainly de- 
termined by those disorder realisations, which move fast. 
In order to determine the scaling exponent /? of v, it is 
thus our task to work out a quantitative description of 
the motion of a particle in a certain disorder realisation 
during a period of time between two collapses of force- 
free points. The temporal distance between two jumps 
of the leftmost force-free position is approximately given 
by 



T = v~\ 



(19) 



because this is the time needed to travel through a cor- 
related region of the disorder (which is of length 1). We 
denote the distance to the new leftmost intersection point 
zi(t) by 



6{t) = z{t) - Zl {t). 



(20) 



Note, that by definition 6(t) is negative. Now, Eq. (fTTj) 
yields the identity 



= (z(t) -vt) = - I dt (zi(t) + 6{t) - vt) 



1 

T 



(21) 



zi (0+) 



Zl (t) 







n 





- zo 

*l(t) 



to 



T 



FIG. 4: Illustration of the motion z(t) in between two jumps 
in the case of smooth disorder. 



several steps. First of all, we note that at t — when Z\ 
and Z2 merge, the relation 



9'(zo) 



(23) 



holds. For t > 0, we can expect that z(t) is still close to 
Zo, so we can expand ((9j) around zq. Writing 



C(t) = z{t) - z 
and using (|23[) , we obtain 



dt( = Tvt 



g"(z ) 



C 2 + o(C 3 ), C(o) = o. 



(24) 



(25) 



For small £, we can neglect the second term on the right- 
hand side and obtain 



C(t) « rut 2 /2. 



(26) 



On time scales £ > to = 2[r«^"(^o)] - ^ the second term 
on the righthand side of ([25]) dominates the time evolu- 
tion and we obtain 



Using Eq. (fl"5j) . we obtain from (|2Tj) 



ft, — h 



1 

T 



dt (0(t)) . 



(22) 



The integral on the righthand side of (|22p depends on 
the velocity. But, in order to use Eq. (f2"2"j) to determine 
the scaling exponent, we have to describe the interface 
motion for t > 0, i.e. in the region of the fast motion 
between the previous and the new force free position. 
We are going to do this separately for the two types of 
disorder considered in this paper. 



1. Disorder with a smooth correlator 

The motion of the interface position after the collapse 
of the two leftmost force-free points is best analysed in 



at) 



tj 



:tn 



g"(z )(t d -ty 2 
Clearly, this result can only be used until a time 

2 



t d 



g"(z y 



(27) 



(28) 



for which C(^i) ^ 1 since we made an expansion in £. It 
shows, however, that for i ^ t < t\ the coordinate z 
increases rapidly until it comes close to the new leftmost 
minimum Zi(t). For t > t\, 9(t) is already close to zero 
and therefore gives only higher order contributions to the 
righthand side of Eq. (j2"2"j) . The motion in between two 
jumps is sketched in Fig. [4] 

Now, we are going to evaluate the integral over 9(t) 
that occurs in Eq. (|22"]) . The equations d20j) and pi)) 
relate 6{t) and £(*) 



6(t)=C(t)+z -z 1 (t). 



(29) 
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The time dependence of z\ (t) can be estimated from Eq. 
IfTT]) as 



d tZl (t) = v + T- 1 g'(z 1 (t))d t z 1 (t) 



Tv 



Tv 



T-g'( Zl (t)) T 



+ 0(v 2 ), 



(30) 



where we have introduced T = [r — g'(zi(0+))] for nota- 
tional convenience. Since zi(0+) is a stable fixed point, 
we have T > 0. Using 



Tv 

z - zi(t) ~ e(o) - —t, 



(31) 



we obtain 

to 



J dt0(t) = dt 



C(t)+<?(0)- y* 



__+tofl(o)- Y - 

=2«-ifl(0)[T fl "(« ! b)]-i+O(l). 
Further, for to < i < ti, we have 



(32) 



ti 

to to 



\n^—^ + e(o)(ti-to)- 



g"{z ) t d -t 
Tv t\ - 4 
Y 2 

=v-h(p)\Tg"(zo)]-v +0(\nv). (33) 

As we have already said, the integral over the remaining 
time interval [fi;T] contributes to 0(1) only. Thus, up 
to orders 0(vhxv), from ([32")) and ([33]) the expression on 
the righthand side of Eq. (|22|) follows as 

T 

- | dt <^(i)> ^ - W §3 (|e(0)|[r 5 "(^ )]-*) • (34) 



From 



and (J3JJ, we obtain therefore 

(h - M 3/2 



(35) 



i.e. /3 = 3/2. 



2. Disorder with a cu. 



correlator 



As we have mentioned before, if A(z) has a cusp singu- 
larity, the typical disorder force realisation exhibits dis- 
continuous jumps, as is depicted in Fig. |3J A moment 
reflection shows, that a merging of two force free solu- 
tions z\ and zi is only possible at such a discontinuity 



T(z-(z))-h 




FIG. 5: Left: This picture corresponds to our assumption for 
case 1, that the new intersection point zi(0+) is left of the 
next discontinuity of the disorder force g(z). It is obvious, 
that the inequality (|36[) has to be fulfilled. Right: A scenario 
contrary to case 1 is possible. However, the basic fact that 
the particle moves from the very beginning with a velocity 
g(z) — [F(z — (2)) — h] = O(l) and therefore needs a time 
to — 0(1) to approach zi(0+), remains unchanged. So does 
the exponent j3. 



of the force field. The requirement, that z\ is a stable 
fixed point entails that such a discontinuity is given by 
an upward jump in the force field. For the calculation we 
have to distinguish several cases. 

Case 1: In this case we assume, that the next stable 
intersection point occurs before the next discontinuity. 
Then, we have the inequality (cf. Fig. [5]) 



T > g'(zo- 



(36) 



It turns out that we have to solve the equation of mo- 
tion in two time regimes. First, close to t = 0, z(t) is in 
the vicinity of zq and we consider again the equation for 



C(t) = z(t) - z . 



(37) 



Now, since the merging of two fixed points occurs at the 
discontinuities of the potential, Eq. (|23j) is not meaning- 
ful, but instead zq fulfills the equation 



r^o = g(zo-) + h. 



(38) 



Using Eq. (f3oj) . it is easy to see that the equation of 
motion for ((t) takes the form 

d t C(t) Ps r(vt - 2 ) + 5(^0+) + (g'(z +) -r)( + h 

= Tvt-(r-g'(z +))((t)+6g. (39) 

Here, 5g = g(zo+) — g(zo—) denotes the jump of g(z) 
which is of order one. Integration of (|39|) gives for short 
times t > 

C(t)^Sgt. (40) 

This result is approximately correct for t < to with 

t = (r -g'(z +))- 1 . (41) 

Note, that due to (|36"|) the time to is always finite and 
positive, in fact generically of order 0(1). For t > to, 
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also the term in Eq. (|39|) proportional to becomes 
relevant. Now, z + ((t ) has to be compared with Zi(0+) 
which is the new leftmost intersection point for t > 0. 
From (JTTJ) we deduce that zi(0+) fulfills the equation 

rz x (0+) « ft + 5 (z +) + ff / (afo+)(*i(0+) - «b), (42) 
from which we obtain 

/i + .g(z +) -g'(z Q +)z 



zi(0- 



C(to) + z a = z(t ). 

(43) 



T-g'(z +) 
[T - g'{z +)}z + Sg 
T-g'(z +) 



In the second step, we have replaced h using Eq. 
Thus, after the time to the particle has reached already 
the new intersection point z\ . 

To determine the exponent /3, we want to employ equa- 
tion ([22"]) again. For t < to, the relevant function 6(t) as 
has been obtained so far reads 



0(t) =z - zi(t) + C(t) nzo- Zl (t) + Sgt. 



(44) 



To approximate the time dependence of z\(t), we expand 
Z\(t) around zi(0+) and get 



z x {t) ~ z x (0+) + zx(0)t. 



(45) 



Here, ii(0) can be deduced from the defining equation 
(HZ}, it follows as ii(0) = Fvti with 



h = (T-g'(z 1 (0))y 



(46) 



Thus, in the regime where 9{t) changes fast, i.e. for t < 
to, we can write 

6{t) Kzo- zi(t) + C(t) ~ * 5 (t - to) - Tvtt t. (47) 

This shows, that for t > t , 6{t) = 0(v). However, the 
time scale to is of the order O(l), and is thus small com- 
pared to T, to *C T. Therefore, it is important to care- 
fully analyse the function 6{t) also for t > t - For t > t n 
we expand around Z\(t) and the approximated equation 
of motion reads 

d t z ~ T{vt -z) + g(zi) + g'(z!)(z -z x ) + h = --U(i), 

ti 

(48) 

where t\ is defined in (|46|) . Then, using (fTT|) and (|20)) . 
we find 



9 t ^~ _ -Q(t)-Tvt x . 
ti 



(49) 



The solution to this equation, matching with equation 
P7)l gives 

6>(*) w -Tvt\ + Tvh {h - t ) e - {t ~ ta)/tl . (50) 
The motion z(t) in between two jumps is sketched in 



zi (OH 



Zl(i) 



Z '(t),'- 



z(t 



- ^0 
Zl (t) 



*0 *0 



FIG. 6: Illustration of the motion z(t) in between two jumps 
in the case of disorder with a cusped correlator. The difference 
between the absence (case 1, solid line) and the presence (case 
2, dashed line) of a discontinuity of the force field g(z) in 
between zo and zi(0+) is the appearance of a sharp kink in 
the curve z'(t) at t = t'j. 



Fig. E 

Now, we can determine (3 using equation (|22[) . In cal- 
culating 



dt (9(t)) 



/Sgt 



Tt\ 



0{v 2 



(51) 



we have decomposed the integral into the intervals ... to 
and to ■ ■ ■ T, respectively. This gives 



h — h 



p 



(52) 



from which we conclude, that in the case of cusped dis- 
order the velocity exponent is j3 = 1. 

Case 2: Now, we have to discuss what can change if 
there is a discontinuity of g{z) in between zq and zi(0+). 
One possible scenario for this case is depicted in the right 
part of Fig. [5j We are going to discuss now, that our 
main result /? = 1 remains unchanged. Indeed, as can 
be concluded from our previous calculation, the essential 
point that lead to the exponent f3 = 1 was the fact, that 
z(t) approaches z\ (0+) on a time scale t which is of order 
0(1). Responsible for this is, that immediately after a 
collapse of the leftmost intersection point, the particle 
starts to move with a velocity of order 8g = 0(1). This 
remains unchanged. In Fig. [6] we have also sketched 
the motion when a discontinuity occurs in between zo 
and zi(0+). The respective quantities in Fig. [S] carry a 
prime. The only effect of the discontinuity that is crossed 
at a time t'j is a singularity of the velocity z'{t) at t = t'j. 
Therefore, the fundamental characteristics of the motion 
remain unchanged. Thus, in case of k jumps of g{z) the 
foregoing calculation remains basically unchanged, apart 
from the fact, that one should now decompose the motion 



in more parts: [OjijJ, [tj 1 ;tj 2 ], 



[t jk ;t ], [h;T\. Of 



course, this consideration changes the prcfactor in Eq. 
()52p . which is, however, anyway beyond our accuracy. 
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0.01 



v(h — h p ), (cusped disorder. num.) o 
fit curve vJh) = 1.32 ■ (ft - h p ) 106 
v(h — h p ). (smooth disorder, num.) + 
fit curve v,(h) = 1.00 ■ (ft - h p ) 151 



0.1 



h — h D 



FIG. 7: The velocity as a function of h — h p for F = 0.67 
in a double logarithmic plot. The numerically determined 
exponent for this measurement is j3 = 1.06±0.08 for cusp-like 
singularity of A (diamonds) and /3 = 1.51 ± 0.08 for smooth 
force correlation (crosses). 




FIG. 8: In the presence of an ac driving field, a velocity hys- 
teresis emerges. In this picture we illustrate these hystereses 
for the cusped and the smooth disorder for T = 0.5, h = 2.0 
and uj — 0.1. The inner hysteresis is traversed clockwise, the 
outer loops are passed through counter-clockwise 



The two exponents (3 = 3/2 for smooth and (3 = 1 for 
cusped disorder are confirmed by our numerical solution 
as depicted in Fig. [7] 



IV. FINITE FREQUENCIES 

In the finite frequency case, the disorder average over 
the solutions to the equation of motion ([7]) forms a hys- 
teresis in the v-h-plane, as is illustrated in Fig. [5] for the 
two types of disorder considered here. The hystereses are 
invariant under the transformation v — > — v and h — > —h. 
This can be explained directly using the equation of mo- 
tion (J7J) and a statistical inversion symmetry. Taking the 
disorder average of ([7]) yields 

























z 









Z r (t b 



z r (t c 



FIG. 9: At h(to) ~ —hp, the particle is close to the rightmost 
force free point Zf(t a )- This intersection point moves, due to 
the change of the zero z r (i) = (2) (t) + h(t)/T. The particle 
is following this point. At a later time t c , when h(t c ) ~ h p 
this force free point has become the leftmost one. 



It is easy to see, that the aforementioned symmetry under 
V — > —v and h — > — h holds true if the probability density 
P(g) (cf. Eq. ©) obeys P(g) = P(g) where g(z) = 
—g(—z). This is obviously the case for our assumption 
of Gaussian disorder (cf. Sec. HI)) . 



A. Qualitative discussion of the motion 

To understand the shape of the hysteresis, we consider 
the motion of a particle for half of a period for the case 
h^> hp and small frequency u ^ T. We start at a time 
t = 0, when h(0) = — h p and the field increases. Then, 
we can expect each particle to be located close to the 
rightmost force free point, i.e. the rightmost solution of 



T(z f (t)-(z) (t))-h(t)=g(z f (t)). 



(54) 



In Fig. [5] it is illustrated, that due to the change of the 
driving field towards larger values, the root of the straight 
line, given by 



(z) (t) + h(t)/T 



moves with a velocity 



v(t) 



r 



(55) 



(56) 



Since h(t) > 0, this velocity is positive although the value 
of the field is still negative. Therefore, also the inter- 
section point Zf(t) to which the particle is connected, 
moves to the right. This fact is observed in the hystere- 
sis loop, illustrated in Fig. [5J Actually, this understand- 
ing allows to estimate the velocity in simple geometrical 
terms. Using the notation explained in Fig. I10| we have 
5g/Szf ~ g'{zf) and thus 



5z f = - 



Sg _ g'(z f )6z f 



from which we conclude 



d t (z) = v = hcosujt + (g(z(t))} 



(53) 



Sz r 



r-g'(zf) 



Sz f . 



(57) 



(58) 
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Zr(t) 



z,.(t + 8t) 
z f (t) z f (t + St) 



right of the leftmost force free position, which happens 
about when h(t) ~ h p , the velocity is negative. Now, the 
same procedure starts in the other direction. 

As a; — > 0, the hystereses approach the depinning curve 
that has been discussed in the previous chapter. In the 
following, we are going to take a closer look on the details 
of this limiting process. 



B. Velocity exponents 



FIG. 10: Illustration for the velocity estimate 



Now, Eq. (J56J) yields 



Sz r 
It 



r-g'(z f )Szf 



r 



St 



r 



Solving the last approximate equality for 5zf/5t, we ob- 
tain 



-h/g'(z f ) 



(60) 



During the motion of Zf(t), other intersection points 
to the left of Zf(t) vanish, and new solutions to the right 
emerge. Finally, when h(t\) ss h p , Zf(t\) has become 
the leftmost intersection point. From approximately this 
time on it happens, that occasionally in some disorder re- 
alisations Zf(t) merges with an unstable fixed point and 
vanishes, so that the particle moves fast in order to catch 
up with the new leftmost force free point. This procedure 
has already been discussed earlier in Sec. IIII1 Since the 
velocity of a particle is given by the difference between 
g(z) and the straight line T(z — (z)) — h, it must fall 
back behind the leftmost intersection point to speed up. 
This can only happen due to the disappearance of force 
free points. Thus, the velocity grows slowly because after 
each jump the particle moves fast and thus approaches 
again the new intersection point. On the other hand, by 
virtue of Eq. (|5B|) . the larger v(t) the faster z r (t) and 
thus also the faster the intersection points move. This 
leads to a positive feedback and entails a strong slope 
when the velocity is large enough such that the particle 
is no longer able to approach a force free point before the 
next jump sets in. Finally, far above h p the particle is de- 
pinned. After the driving force has reached its maximum 
it decreases. Note that the root of the straight line z r has 
now a velocity smaller than v(t), because h is negative. 
Therefore, the particle position z(t) approaches z r (t) and 
slows down. Hence, h is a measure also for the decrease 
of v(t). On approaching h p from above, all particles are 
still dcpinned and hence far enough behind the leftmost 
intersection point, so that the latter has only little in- 
fluence on the motion of the particle and the velocity 
decays with the same slope all the time. Only when v(t) 
has passed below h/T, z r moves in the negative direction 
and thus the intersection points as well. This means, 
that the leftmost intersection point approaches the par- 
ticle before it is pinned. After the particle is a little to the 



First, we want to work out, how 
v ho = \v(h = 0)\ 



(61) 



approaches zero as u — > 0. As we have explained in Sec. 
IIV A[ the particle in each disorder realisation stays close 
to a force- free point, that we have agreed to label zf(t). 
The velocity dtz of the particle is now determined solely 
by the velocity of the force free position Zf that we are 
now going to calculate in a more accurate way than our 
estimate from Eq. ([SU]). Let to be the point in time, at 
which h(to) = 0. On time scales that are small compared 
to the period 0J , we can linearly expand the driving 
field around to 



h(t) ~ -hu(t - t ). 



(62) 



Further, we want to expand (|54[) around zj(to)- For small 
distances in time we can neglect possible changes in the 
velocity and write 



z f(t) -Zf(to)+Vf(t~t ), 



(63) 



where Vf — dtZf is a shorthand notation. Using (j()2"j) as 
well as (z) (t) ~ (z) (to) — V] l() (t — to), we have 

=T[z f (t ) + v f (t - t ) - (z) (to) + v ho (t - t )] + 
hw(t - t Q ) - g(z f (t Q )) - g'(z f (to))v f (t - t ) (64) 
= (t - to) [T(v f + v ho ) +hu- v f g'(z f (to))} (65) 

Since this should hold for small but finite |t — to|, the ex- 
pression in the rectangular brackets has to vanish. Solv- 
ing (|65j) for Vf , taking the disorder average and using the 
self-consistency condition (vf) = — Vh finally yields 



Vh, 



T-([T-g'(z f (to))}-'Y 



(66) 



Since g'(zf(to)) < 0, which expresses the reasonable 
assumption that Zf is a stable force free position, Vh is 
indeed positive, which must be the case by its definition 
(|6T|) . Note, that our derivation so far does not make any 
assumption about the disorder correlator, whence it holds 
for cusped as well as for smooth disorder. In conclusion, 
for lu — > the width of the hysteresis at h = behaves 
as Vh ~ w K with k = 1 for either kind of disorder. This 
exponent is verified by our numerical analysis, cf. figure 

HU 
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FIG. 11: The velocity Vh as a function of frequency. The 
plotted data correspond to numerical measurements at F = 
0.33 for cusped and F = 0.2 for smooth disorder. For both 
types of disorder, the exponents are close to 1 (k c = 0.97 ± 
0.07 and k s = 0.95 ± 0.04) in agreement with our analytical 
derivation given in the main text. 



case in 4 — e dimensions, the exponent fi has been found 
as n = f3/(vz), where v denotes the correlation length 
exponent and z the dynamical exponent [2l| . 



C. The area of the hysteresis loop 

Next, we want to investigate the limiting behaviour 
of the hysteresis area. The physical meaning of the loop 
area can be concluded from the energy balance of an over- 
damped system. For the change of the disorder averaged 
potential energy in time, we find (cf. App. |E| 



8 t (E) = h(t)v(t) - {{d t zf) 



(68) 



Here, h(t)v(t) measures the energy gain through the work 
per unit time that is done by the external field and 
((<9tz) 2 ) measures the energy loss per unit time due to 
dissipation. The area of the hysteresis loop is determined 



o.i 
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Vk p (uj), (cusped disorder, num.) 

fit curve v c (uj) = 0.67 • u 0AS 
Vh p {^S), (smooth disorder, num.) 
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FIG. 12: The velocity Vh p as a function of frequency. The 
plotted data correspond to numerical measurements at V = 
0.33 for cusped and V — 0.2 for smooth disorder. For both 
types of disorder, the exponents are close to 1/2 (fj, c = 0.48 ± 
0.02 and fi s = 0.51 ±0.01). 



Another interesting quantity to look at is 
v h = \v(h = h p )\, 



(67) 



of which we want to work out the limiting behaviour for 
u> — > 0. As h(t) increases further from towards h p , more 
and more of the force free points zj, which the true po- 
sitions in the disorder realisations arc following, become 
the rightmost ones, so that occasionally jumps occur. On 
closely approaching h(t) = h p the dominant contribution 
originates from these jumps, which severely affects the 
exponent, so that Vh p ~ ^ with fj, ~ 1/2, as can be 
inferred from our numerical analysis, shown in figure [121 
This exponent is again independent of the shape of the 
disorder correlator at the origin (smooth or cusped). An 
analytical derivation of this exponent is much more com- 
plicated than it was the case for k and in fact we did 
not find a rigorous prediction. For the finite dimensional 



^hyst(w) = fv(t)dh= / v(t)h(t)dt 



(69) 



This means, the loop area denotes the integrated change 
in work per unit time due to the change of the external 
field. 

Note, that in the case of a double hysteresis (which 
occurs for large h, when the motion of the system over 
one period extends on average over more than one valley 
of the disorder potential), the area is given by the area 
of the inner hysteresis minus the area of the two outer 
hystereses (cf. Fig. [8]). Formally, this is because the inner 
hysteresis is traversed clockwise, whereas the outer loops 
are passed through counter-clockwise. Physically, this 
can be understood as follows. Starting from h = at the 
branch of increasing h(t), the external field works against 
the potential gradient due to elastic energy and disorder. 
On going over into the regime of the outer loops, sliding 
behaviour sets in and thus the potential energy, stored so 
far in the system, adds to the work done by the external 
field. This fact is responsible for the steep slope at the 
beginning of the outer loop. In other words, during the 
period in the outer loops, the external field does not any 
more work against a potential gradient, but together with 
the potential energy the system is accelerated. 

To work out the hysteresis loop area as u — > 0, we 
distinguish three cases. 

(a) h <C hp. In this situation, the hysteresis con- 
sists of a single loop. The outer loops, visible in Fig. 
[51 are absent. We expect the loop area to be given by 
^4-hyst ~ Vh h ~ uj k . Indeed, our numerical solution shows 
that the of the hysteresis vanishes proportional to w, 
independent of the type of disorder correlator, as shown 
in Fig. [HI 

(b) h = hp. For this case, the hysteresis loop is still 
single (no double hysteresis) and the hysteresis area de- 
creases with the frequency as ^4hyst ~ w 82 (cf. Fig. IT4"|) . 
still independent of the disorder correlator. 
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FIG. 13: The area of the hysteresis loop is plotted as a func- 
tion of the driving frequency in case h <C h v . For small uj the 
area is diminished proportional to uj (exponent 0.98 ± 0.01), 
as expected. 
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(c) h ^> h p . Now, we face the situation of a double 
hysteresis and moreover, the behaviour of the hysteresis 
area as u —> now depends on the shape of the disorder 
correlator. We find Ah ys t ~ w Q with a ~ 0.67 w 2/3 for 
cusped and a ~ 0.75 = 3/4 for smooth disorder. This is 
shown in figure 1151 

So far, our results suggest that the scaling exponents 
are insensitive to the nature of the disorder correlator as 
long as the force amplitude does not exceed the threshold 
hp. These findings seem to milden the non- universality 
conclusion by Fisher ||, who considered the response of 
a charge density wave system to an ac force in addition to 
dc driving. He distinguished different distributions of the 
random amplitude (pinning strength) of the disorder po- 
tential in addition to a random phase, and found a strong 
dependence of the behaviour on the type of disorder both 
above and below threshold. 

For large frequencies, the area of the hysteresis loop 
vanishes as well. Above a certain crossover frequency, 
which depends on h, the motion of the particle is re- 
stricted to one minimum of the potential. Thus, for large 
enough frequencies we can approximate the potential by 
a harmonic one, such that the equation of motion for the 
disorder averaged position becomes 



v(t) = h cosujt — w z(t), 



which has the solution 
wujh 



v 



(*) 



W 2 + UJ 2 



— sin H cos cut 

w 



(70) 



(71) 



o.oi 



0.1 



Thus, using Eq. (|69l) we find for the hysteresis loop area 



FIG. 14: The area of the hysteresis loop is plotted as a func- 
tion of the driving frequency in case h ~ h p . For small uj the 
area vanishes with an exponent 0.82 ± 0.01, independent of 
the type of disorder. 



^hystH = J V(t)h(t)dt 




1" 

;— ~ UJ 



(72) 



where the last expression gives the asymptotics for large 
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FIG. 15: The area of the hysteresis loop is plotted as a func- 
tion of the driving frequency in case h~3> h p . The diminution 
of the area with cj can be described as a power law with dif- 
ferent exponents for cusped (a = 0.67 ± 0.03) and smooth 
disorder (a = 0.75 ± 0.04). 



The decay of the hysteresis loop area for large and 
small frequencies of the driving force requires the exis- 
tence of a maximum. This maximum is found to be pro- 
portional to the resonance frequency of the typical dis- 
order potential wells uj r = juq/1, which equals 1 in our 
units. The proportionality factor is of order unity, and is 
found different for small driving fields (single hysteresis) 
and large drivings (double hysteresis). 



V. CONCLUSIONS 

For the mean-field theory of driven elastic manifolds 
in disordered systems, we have worked out the scaling 
behaviour of the velocity as a function of the dc-driving 
force close to depinning by extending Fisher's [1] argu- 
ments for charge density waves. The scaling exponents 
are found to be different for disorder with smooth and 
cusped correlator. Our analytical findings arc supported 
by a numerical treatment. 
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Furthermore, we have investigated the small-frequency 
behaviour of quantities that characterise the velocity hys- 
teresis in case the system is exposed to an ac-driving. We 
found that the frequency scaling exponents of the rema- 
nent velocity Vh and the velocity at the depinning field 
Vh p do not depend on the presence or absence of a cusp- 
singularity at the origin of the disorder correlator. This 
also holds for the frequency exponents of the hyteresis 
loop area as long as the amplitude of the driving does 
not exceed the depinning force h p . For force amplitudes 
above the depinning threshold, our numerical treatment 
yields different exponents for smooth and cusped corre- 
lators. 
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Appendix A: Applicability of the mean field theory 

Mean field theory is valid above the upper critical in- 
terface dimension D c , where the interface is not rough 
and there is no depinning transition for weak disorder. 
Below this critical dimensionality, the interface is soft 
enough to adapt to the random potential. Hence, it be- 
comes rough and is pinned for small driving forces. An 
FRG analysis of the depinning transition shows, that the 
exponents become non-classical and the disorder corre- 
lator develops a cusp singularity [13] • The exponents 
can be expressed in terms of an e-expansion, where 
e = D c — D, and the zeroth order O(e ) is given by the 
critical exponents of mean-field theory with a cusped dis- 
order correlator [3~J ]. 

For systems with short range interaction it turns out 
that D c = 4, so interface dimensionalities above D c seem 
to be a purely academic problem. There are, however, 
systems with long range forces for which the upper criti- 
cal dimension D c is reduced. For example, for interfaces 
in systems with long range dipolar interaction it has been 
argued in Ref. [33| that the upper critical dimension de- 
creases to D c = 3. A further decrease of the critical 
dimension for the statics of interfaces is achieved for do- 
main walls in ferroelectric materials with a piezo effect in 
the paraphase. An example for such a material is KDP. 
A peculiarity of such systems is that interfaces between 
different ferroelectric domains are allowed to be oriented 
along some distinguished orientations in the crystal only. 
This is because ferroelectric domain walls are at the same 
time fcrroclastic domain walls which have to fulfill certain 
mechanical compatibility relations [HI, Looking at 
the statics, one finds that in the presence of random field 
disorder the interface is not rough for D > 2 [H, [13] ■ In 
the following, we are going to take a look on the impact 
of these long range forces on the dynamical properties. 



Firstly, we note that our model equation of motion ([T]) 
can be written in a more general fashion as 



—d t z(x, t) ~u a g(x, z) + h(t) 

7 



where 



+ J d d x'T cl (x-x')z(x',t), 

r cl (x) = / e lfex r cl (fc). 



(Al) 



(A2) 



In (TTJ we assumed that only short range forces are present 
and hence 



r c i(fc) = r fc 2 



(A3) 



For the systems with long range forces, discussed above, 
the functional form of r c i(fc) is different from (|A3j) . In 
systems with dipolar interactions, one finds [|[ 



r c i(fc) = r fc 2 + a 



d-2 



[(Ho) 



2-d 



1 



(A4) 



where A measures the strength of the dipolar interaction 
and £o 1S of the order of the lattice spacing. Here, k a is 
the component of k along the orientation of the dipoles. 
A further increase of the domain wall stiffness is achieved 
for domain walls in ferroelectric materials with a piezo ef- 
fect in the paraphase. Integrating out the elastic degrees 
of freedom in such a (3d) system leads to an interface 
Hamiltonian of the form 13611 



r el (fc) = r fc 2 + \f 



ck. 



(A5) 



Here, c is a constant that depends on the strength of the 
piezo coupling as well as on the elastic constants and A 
depends on the strength of the dipolar coupling, the piezo 
coupling coefficient and the elastic constants. Moreover, 
k a denotes the component of k in the direction of the 
polarisation. 

In the following, we focus on the case of a dc-driving 
field h(t) = h. For the perturbative determination of 
the velocity, it turns out to be useful to go over to a 
co-moving frame z(x,t) = vt + £(x, t). Expanding the 
correction t) in a power series in uq, 



C(x,t) = J2 U oCn(x,t) 1 



(A6) 



one can solve Eq. (|A1[) pcrturbativcly by comparing 
equal orders of uq. Following the computation in Ref. 
[l8| (which is similar to that presented in App. [C] be- 
low), we obtain a self-consistent perturbative expression 
for v (th e following equation corresponds to Eq. (29) in 
Ref. p||) 



v — 7/1 



q Jk 



; ''</ 



r e i(fc) 

iqA q ul(T el (k)+i^) 



q J k 



[r c i(fc)P 



(A7) 
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FIG. 16: An example for a randomly generated disorder re- 
alisation with the method that leads to a statistics with a 
smooth correlator. Besides the force g(z), also the disorder 
potential, determined from g(z) = —V'(z) is visualised. The 
potential has a smooth shape. 
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FIG. 17: An example for a randomly generated disorder re- 
alisation with the method that leads to a statistics with a 
cusped correlator. Besides the force g(z), also the disorder 
potential, determined from g(z) — —V'(z) is visualised. It 
can be seen that most of the extrema of the potential reveal 
a cusp singularity. 



Here, A q denotes the Fourier transform of the disorder 
correlator defined in Eq. ©. For v — > 0, our choice (|A3|) 
for r e i(fe) leads to a divergence of the integral in (|A7|) 
when D < 4. The divergence stems from the contribution 
of small A;. In contrast, for a model the elastic interaction 
of which assumes the form given by Eq. (|A5[) , the integral 
in Eq. (|A7|) is convergent as long as D > 2 and exhibits 
only a logarithmic divergence in D = 2. Thus, in such a 
system mean field results may apply for two-dimensional 
domain walls with small logarithmic corrections. 
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Appendix B: Numerics 



-0.2 



In this appendix, the details of our numerical methods 
shall be explained briefly. Numerically, we have solved 
the equation of motion ([7]) for 20000 different disorder 
configurations using the well established classical Eulcr 
method. Before embarking on the details of the proce- 
dure in the two cases of dc and ac driving, respectively, 
we want to explain how the different realisations of the 
disorder forces are achieved. 

The random force associated to the smooth disorder 
correlator is generated by concatenating straight lines at 
distance 1. The values of the disorder force at the con- 
catenation points are randomly chosen out of a bounded 
interval [— m; m]. The bounds to for this interval are de- 
termined such that the correlator fulfills A(0) = 1. Wc 
found to ~ 2.1. Moreover, the position s of the con- 
catenation point closest to the origin z = has been 
determined randomly out of the interval \—\\ |]. This 
method has already been used earlier by A. Glatz [38| . 
A sample configuration is depicted in Fig. [16J 

To realise disorder forces with Gaussian statistics that 
arc correlated with a cusp singularity at the origin, we 
have also used straight lines, the extension in z-direction 
being again 1. In contrast to the former case, these lines 



FIG. 18: The numerically determined correlators for the dis- 
order forces generated by the two different methods described 
in the text. Either disorder correlator decays on a length scale 
of order unity and obeys A(0) = 1. 



are now concatenated discontinuously. The values of g(z) 
at both endpoints of a segment are determined randomly 
from a bounded interval [—to; to]. The size of this interval 
is again determined in such a way that the disorder cor- 
relator obeys A(0) = 1 with the same result of m ~ 2.1 
as before. Also the position s € [— |; 4] of the jump clos- 
est to z = is created randomly. An example for such a 
configuration is shown in Fig. 1171 In App. |D]wc discuss 
the latter generation technique in more detail and also 
derive the force correlator. 

Moreover, we have verified numerically, that the as- 
sumed shapes for the disorder correlators are reflected 
by our two generation techniques. The result is shown in 
Fig. [THJ As required, one correlator is smooth and the 
other shows a cusp singularity at the origin, perfectly in 
agreement with the analytic result (cf. Eq. (|D8|) below). 
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In both cases, the correlations decay to zero on a length 
scale I ~ 1. 

To solve the equation of motion with a constant driving 
force for several values of T and h, we have chosen the 
initial condition v(0) = and z(0) = 0. The time steps 
are chosen to be of size 2 -10 and we have simulated the 
equation of motion for 2 14 such time steps. To avoid that 
our results arc spoiled by transience effects, we have only 
taken into account the values of v for the last 2 10 time 
steps. 

Our results that we have obtained for the ac driving 
case rely on the initial conditions z(0) = 0, v(0) = h 
with h(t) = h cosujt. Before making any measurements 
we have been waiting for at least 2 periods for transience 
effects to diminish. The values for v(h = 0) and v(h = 
h p ) are measured then over 3 periods with 2 datapoints 
per period. 

In all cases where we have numerically determined im- 
portant exponents, we have given error estimates. These 
include the statistical deviations. 

Appendix C: Perturbation theory 

To set up the perturbation expansion for our equation 
of motion ([8]) 

d t z = T({z)-z) + h + g(z), (CI) 

we write z(t) = (z(t)) + dtz = v + dt(, so that 
(£) = 0. This yields two coupled differential equations 

v = h+(g((z)+0) (C2) 
(d t +rK = g((z)+0-(g((z)+C)). (C3) 

The differential operator on the left hand side of (|C3[) 
has the fundamental solution 

{d t +T)G(t)=S(t), (C4) 

G(t) = Q(t)e~ rt . (C5) 

Now, for large enough T, we can expect ((t) <C 1 and 
perform a Taylor expansion of the disorder force 

^> + o = E^c". (06) 

Thus, (|C2j gives us 

v = h + (g((z))+g>((z))C + ...) 
= h + (g'((z))C(t) + ...) (C7) 

For ((t), we use the lowest order result from the iteration 
([C3)) combined with ([C6} 

OO 

C(t) = J dt' G(t - [g((z(t')) + C) - (g((z(t>)) + C)>] 
/ 

= J dt'e- T ^\g{{z{t'))) + ...}, (C8) 

to 



where t$ is the time at which we fix our initial value 
problem. Sending to ~ ^ — oo, we obtain 

v = h+ J dt'e-^- t 'Hg\{z{t)))g{{z{t l )))) + ... 

— OO 

/ 

= h+ J dt'e- r ^A'[(z(t))-(z(t'))}. (C9) 

— oo 

A change of the integration variable t' — > t—t' and noting, 
that (z(t)) = vt then immediately yields the result given 
by CGI. 

Appendix D: Generation of disorder forces with a 
cusped correlator 

In this appendix, we propose a generation technique for 
disorder forces with a cusp singularity in the correlator 
and prove, that the cusp is indeed present. The method 
described here has also been employed in our numerical 
analysis, cf. App. |B] 

Consider a function g(z), that is constructed as follows. 
We decompose the z-axis in intervals Ii — [z^z^i] of 
length 1. The starting point of each interval Ii is given 
by 

Zi = i + s, (Dl) 

where i is an integer (the label for the interval) and s is 
a random number, uniformly distributed in the interval 
[— i; i]. The function g(z) shall now be given piecewise 
for each interval Ii as a straight line. This line is deter- 
mined by the left boundary point (zi,ai) and the right 
boundary point (zj+i, (3i), where en and are chosen 
randomly out of a bounded interval [— m; m}. The value 
of m will be determined further down in such a way, that 
A(0) = 1. 

In more mathematical terms, a specific realisation of 
the function g(z) is specified by the random number 
s G [—5 : 5] an d two stets of boundary values {cti} and 
{/3i}. Using the indicator function K(z; Zi, Zj+i) for each 
interval Ii, defined by 

*(*5*.*h) = {; Zi ol Z erwi Z se +1 ' ( D2 ) 
the function g(z) is explicitly given by 

g(z, [s; {ai};{f3i}}) = ^ K(z; z h z l+1 ) 

i 

x[(Pi-ai)(z-i-8) + ai\. (D3) 

Examples for typical configurations of g(z) are depicted 
in the Figs. [3land[T71 An illustration for a single segment 
is sketched in Fig. Q1JJ In the following, we are going to 
show, that disorder forces g(z) given by (|D3[) fulfill our 
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i + s 



i + l + s 



FIG. 19: A segment of the disorder force with a cusp singu- 
larity in the correlator. The values a,i and /3j are chosen out 
of a bounded interval. 



requirements (g(z)) = and (g(z)g(z')} = A(z — z'), 
where A(z — z 1 ) obeys A(0) = 1, decays to zero over a 
length scale of order 1 and shows a cusp singularity at 
the origin. Straightforwardly, we find 



ds 



d/3da 
Am 2 



[(/3 - a)(z - i - s) + a] = 0. 



(D4) 



To calculate the second moment (g(z)g(z')), we have to 
distinguish two cases. If \z — z'\ > 1, the points z and z 1 
have to belong to two different intervals z G I k , z 1 G I k > , 
because their distance is then larger than the size of an 
interval and hence 

9(z)g(z') = ^2 K ( z > z *' [(A - oii)(z - Zi) + ati] 
i,j 

x K(z';Zj,z j+ x)[(Pj - otj){z' - Zj) + atj] 
= [(/3k - a k )(z - k - s) + a k ] 

x[(J3 k >-a k >)(z'-k'-8) + a k ,]. (D5) 

This gives for the correlator 



(g(z)g(z')) = / ds 



d(3 k df3k>da k da k > 
16to 4 

_ I — m 

2 

x [(J3 k - a k )(z - z k ) + a k ] 
x [{j3 k > - a k >)(z' - z k >) + Oik>] 
=0. (D6) 

On the other hand, in case \z — z'\ < 1, those realisations 
which do not have a jump in between z and z', i.e. for 
which z and z' belong to the same interval I k , give a 



finite contribution. These correspond to values of s, for 
which 



£ = \z — z | + min(z, z ) — i — s 



(D7) 



obeys \z — z'\ < ( < 1. Instead of integrating over s, it 
is easier to integrate over £ which gives us 



x [(fa - a k ){C - \z - z'\) + a k ] 



2m 1 



l-°-\z-z'\+ L -\z-zf 



(D8) 



As required, the correlator exhibits a cusp singularity at 
the origin and decays to on a length scale I = 1. To 
fulfill A(0) = 1, we have to take m = 3/\/2 ~ 2.1, which 
is in agreement with our numerical technique, described 
in App. [Bj 



Appendix E: Derivation of the potential energy 
balance 



The equation of motion can be considered to follow 
from a Hamiltonian 



d t z 



6H[z] 
~5z~' 



with 



H[z] = E[z] - h{t)z{t) 

E{z] = ^(z(t)-(z(t))) 2 + V(z). 



(El) 

(E2) 
(E3) 



Here, E[z] denotes the total potential energy, V{z) is the 
disorder potential related to g(z) via g(z) = —V'(z). The 
change of the potential energy in time thus follows as 

d t E[z] = T(z - (z))(d t z - v(t)) - g(z)d t z 

= T((z) - z)v(t) - [T((z) -z)+ g(z)] d t z. (E4) 

Using Eq. (JT]), we can replace the term in the square 
brackets and obtain 

d t E[z] = - z)v(t) + [h(t) - d t z] d t z. (E5) 

Taking the disorder average readily yields the result, 
stated in Eq. (|68]) . 
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